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ABSTRACT 

The gravitational lensing properties of cosmological halos depend upon the mass distri- 
bution within each halo. The description of halos as nonsingular, truncated isothermal 
spheres, a particular solution of the isothermal Lane-Emden equation (suitably modi- 
fied for A ^ 0) , has proved to be a useful approximation for the halos which form from 
realistic initial conditions in a CDM universe. The nonsingular TIS model reproduces 
many of the quantitative features of the N-body results for CDM halos, except in the 
very center, where CDM N-body halos show density profiles which vary as p ~ r _Q , 
a > 1, instead of a small flat core. Possible discrepancies between these cuspy halo 
predictions of the CDM N-body simulations and observations of the inner mass pro- 
files of dwarf and LSB disk galaxies based upon their rotation curves and of clusters 
based upon strong lensing measurements have led to a search for other diagnostics. 
A description of the lensing by TIS halos would be useful in this regard, as a self- 
consistent model for CDM halos in a proper cosmological context, nonsingular but 
otherwise consistent with the CDM N-body results. 

We derive here the basic lensing properties of individual TIS halos. For compari- 
son, we also consider three singular profiles: the Navarro-Frenk- White density profile, 
the singular isothermal sphere, and the Schwarzschild lens. For all profiles, we compute 
the critical curves and caustics, the image separation, the magnification and bright- 
ness ratio, the shear, the time delay, and the average shear inside the tidal radius. 
This provides tools for studying the statistical properties of lensing by TIS and other 
lenses in the context of a theory of cosmological structure formation. 

Key words: cosmology: theory - dark matter - galaxies: clusters: general - galaxies: 
formation - galaxies: halos - gravitational lensing 



1 INTRODUCTION 

The gravitational lensing of distant sources has in recent 
years become one of the most pow erful tools in observa- 
tional cosmolog y (see, for example. | So ucail||200ll and ref- 
erences therein; iBartelmann fc Schneiderl l200ll and refer- 
ences therein). Since the effects of gravitational lensing de- 
pend upon the redshift of the source, the cosmological back- 
ground, and the distribution of matter in the universe, they 
can be used to constrain the cosmological parameters and 
the primordial power spectrum of density fluctuations from 
which structure originates. In addition, many of the effects 
produced by gravitational lenses, such as image multiplic- 
ity, separations, and time delay, depend strongly upon the 
matter distribution inside the lenses. Hence, measurements 
of these effects can provide a unique tool for probing the 
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matter distribution inside collapsed objects like galaxies and 
clusters, providing the only direct measurement of their dark 
matter content, and constraining the theory of their forma- 
tion and evolution. 

Until recently, the internal structure of halos adopted 
in lensing studies was generally some gravitational equi- 
librium distribution, either singular or nonsingular (e.g., 
King model, singular isothermal sphere, pseudo-isothermal 
sphere), not necessarily motivated directly b y the theory 
of cosmological halo formation (see, e.g.. lYoung et alJ 



1980; Turner. Ostr iker. fc Gotd 11984 Iffins haw fc Krauss 
198ft iNaravan fc White! ll988MBlandford et al l I l991 
^g^^^^^^ro^S 1^921 |Ko^hanek| ^95 



Premadi. Martel. fc Matznejll998l : iPremadi et alJl20oTglt 
Paisin^^M^n^Olll . As the theory of halo formation in 
the CDM model has advanced in recent years, however, 
the halo mass profiles adopted for lensing models have 
been refined to reflect this theory. Numerical simulations 
of large-scale structure formation in Cold Dark Matter 
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(CDM) universes predict that galaxies and clusters have 
a singular density profile which approaches a power law 
p oc r _n at the center, w i th the exponent n ranging from 
1 to 1.5 llCole fc Lacevl Il996t iNavarro. Frenk. fe White! 
19961. fl997l: iTormen^ Bouchet. fc White! Il997 

d_Li 



Fukushiec & Makin ol Il99l koOlaUbj: iMoore et al 
19991: iHuss Jain fc Steinmetzl I l999t iGhigna et all 



Jing fc Sutoll2000l: iKlypin et alJ l200a iPower et alJ 112002ft 

These results are in apparent conflict with observations of 
rotation curves of dark-matter-dominated dwarf galaxies 
and low surface brightness gala xies , which favor a flat- 
density core fcf. iPrimack et aL | |l999t iBurkert fc SiHdfl999l: 
IMoore et aT]ll999l: lMoorell200ll) . On the scale of clusters 
of galaxies, observations of strong gravitational lensing of 
background galaxies by foreground clusters also favor the 
presence of a finite-density core in the centers o f clusters 
(see, e.g.. lTvson. Kochanski fc DeirAntoniolll99Sl) . 

Several possible explanations have been suggested in 
order to explain this discrepancy. The rotation curve data 
might lack sufficient spatial resolution near the center to 
distinguish unambiguously between a density profile with 
a flat-density core and one w ith a singular profile (e.g. 
van den Bosch fc SwaterslboOll) . Attempts have also been 
made to improve the numerical resolving power of the 
simulations to obtain a more accurate determination of 
the slope of the predicted density profil es at small radii 
fe.g. lMoore et alJll999t iPower et al"ll2002l) . However, if the 
flat-core interpretation of the observations and the singular 
cusps predicted by the numerical simulations are both 
correct, then the simulation algorithms may be ignoring 
some physical process which would, if included, serve to 
flatten the halo density profiles at small radii relative to 
the results for purely gravitational, N-body dynamics of 
cold, collisionless dark matter, while retaining the more 
successful aspects of the CDM model. For example, gasdy- 
namical processes which involve dissipation and radiative 
cooling, and perhaps energy-re lease feedback associated 
with st ars and quasars (see, e.g. iGelato fc Sqmmer-Larsenl 
119991: lEl-Zant. Shlosman. fc Hoffman! l200ll) or mod- 
ifications of the microscopic properties of the 
dark matter, such as self-in te racting dark mat- 
ter iSpergel fc Steinhardl |2000l: lYoshida et all l2000t 
lHannestad. fc Scherrerl 1200(1 iDave et alJ l200ll) warm 
dark matter ~^^^^^^^^ Re^s^^^Valerizuelal 1200' " 



Sojnmer-Larsen 'fe^DoIgovl ^^^)' , fluid dark matt er JPeeblea 
2000), fuzzy dark mat ter l|Hu. Barkana. fc Gruzinov| [2000'). 
decaying dark matter j Ceg^o" ^. an nihilating dark matter 
iKapling hat 1 Knox 1 fc T urneri ~ l2000l) . or repulsive dark 
matter jGoodmanll200Cl) . all have the potential to lower the 
central density of halos and possibly reconcile simulations 
with observations. 

Lensing by the two kinds of halo mass profiles, sin- 
gular versus flat-core, will be different. This has led to 
attempts to predict the differences expected if the halos 
have the singular cusp of the Navarro-Frenk- White (NFW) 
or Moore profi l es or else a profile with a flat core (e.g . 
Kochaneklll995l: iKeeton fc Madaull200lt iRusin fc Ma| |2001: 



Wvithe. Turner. fc Sperge]||200ltTTakahashi fc ChibalboOl : 
Li fc Qstrikerll2002l) . Singular profiles like that of NFW are 
physically motivated by the N-body simulations, and the 
latter have been used to place these halo profiles empirically 
in a proper cosmological context which permits statistical 



predictions for the CDM model. The nonsingular profiles 
which have been adopted to contrast with these singular 
ones, however, are generally no more than parameterized, 
mathematical fitting formulae, with no particular physical 
model to motivate them or put them in a proper cosmolog- 
ical context. 

We have developed an analytical model for the post- 
collapse equilibrium structure of virialized objects that 
condense out of a cosmological background universe, ei- 
ther m atter-dominated or fla t with a cosmological con- 
stant l[ Shapiro. Iliev. fc Ragal Il999l hereafter Paper I; 
llliev fc Shapiroll2001al hereafter Paper II). This Truncated 
Isothermal Sphere, or TIS, model assumes that cosmologi- 
cal halos form from the collapse and virialization of "top- 
hat" density perturbations and are spherical, isotropic, and 
isothermal. This leads to a unique, nonsingular TIS, a par- 
ticular solution of the Lane-Emden equation (suitably mod- 
ified when A / 0). The size r± and velocity dispersion ov 
are unique functions of the mass M and formation redshift 
z C oii of the object for a given background universe. The TIS 
density profile flattens to a constant central value, po, which 
is roughly proportional to the critical density of the universe 
at the epoch of collapse, with a small core radius ro ~ rt/30 
(where a\r = AnGparQ and rp = rKing/3, for the "King ra- 
dius" r-Ring, defined bv lBinnev fc Tremaineiri987l. p. 228). 

Even though the TIS model does not produce the cen- 
tral cusp in the density profile of halos predicted by numeri- 
cal CDM simulations at very small radii, it does reproduces 
many of the average properties of these halos quite well, 
suggesting that it is a useful approximation for the halos 
w hich result from more re alistic initial conditions (Papers I, 
II: llliev fc Shapiroll2001bl and references therein) . In partic- 
ular, the TIS mass profile agrees well with the fit by NFW 
to N-body simulations (i.e. fractional deviation of ~ 20% or 
less) at all radii outside of a few TIS core radii (i.e. outside 
a King radius or so). It also predicts the internal structure 
of X-ray clusters found by N-body and gasdynamical simu- 
lations of cluster formation in the CDM model. For exam- 
ple, the TIS model reproduces to great accuracy the mass- 
temperature and radius-temperature virial relations and in- 
tegrated mass profiles derive d empirically from the sim- 
ulatio ns of cluster formation fevrard. Metzler. fc Navarrol 
1996). The TIS model also successfully reproduces to high 
precision the mass-velocity dispersion relation for clusters 
in CDM sim ulations of the Hu bble volume by the Virgo 
Consortium jEvrard et al.ll2C)02l) . including its dependence 
on redshift for different background cosmologies. The TIS 
model also correctly predicts the average value of the virial 
ratio in N-body simulations of halo formation in CDM. 

The TIS profile matches the observed mass profiles of 
dark-matter-dominated dwarf galaxies. The observed rota- 
tion curves of dwarf galaxies are generally well fit b y a den- 
sity p rofile with a finite density core suggested bv feurkerd 
(1995), given by 



p(r) 



Po,b 



(r/r c + l)(r 2 /r 2 c + l) 



(1) 



The TIS model gives a nearly perfect fit to this profile, with 
best fit parameters Po,b/po,tis = 1.216, r c /ro,Tis = 3.134, 
correctly predicting the maximum rotation velocity u ma x 
and the radius r max at which it occurs. The TIS model 
can also explain the mass profile with a flat density core 
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measured by ITvson. Kochanski fc DeirAntoniol jl998l) for 
cluster CL 0024+1654 at z = 0.39, using the strong grav- 
itational lensing of background gal axies by the cluster t o 
infer the cluster mass distribution dShapiro fc IlievI l200Cft . 
The TIS model not only provides a good fit to the projected 
surface mass density distribution of this cluster within the 
arcs, but also predicts the overall mass, and a cluster veloc- 
ity dispersion in cl ose agreement with t he value ov = 1150 
km/s measured bv lDressler et al.1 lll999D . 

Therefore, the TIS model can be applied to clusters of 
galaxies or dark-matter-dominated dwarf galaxies, for which 
baryonic processes have not significantly modified the mass 
profile. At the intermediate scale of large galaxies, where the 
central density profiles may be baryon-dominated, following 
radiative cooling by the baryon gas, the central profiles may 
differ from that of the halo of dark matter and baryons which 
would otherwise form in the absence of radiative losses by 
the baryonic gas. The TIS model ignores such processes, just 
as the empirical density profiles (NFW, Moore, . . .) based 
on gravity-only N-body simulations do. Hence, these profiles 
might not be directly applicable to galaxy-scale objects for 
which the mass profiles have been significantly affected by 
such baryonic processes. 

Several authors have stu died the effect of lensing by ha- 
los w i th a flat-density core ljaroszvriskilll99ll: IJaroszvnskil 
I l992l: iKochanekl 11993: IPremadi. Martel. fc Matznerl Il99a 
IPremadi et alJl2001al a> or by NFW or Moore profiles that 
have been generalized, so that the inner slope of the den- 
sity profile is arbitrary |K eeton fc Madau|200l|:|Rusin fc Mai 
l200ll:IWvithe. Turner, fc Spergel200ll:lLi fc OstrikeJ2002l) . 
These particular density profiles are essentially mathemat- 
ical conveniences without physical motivation. There is no 
underlying theoretical model in these cases that was used to 
predict the value of the core radius or the departure of the 
inner slope of the density profile from the value found by 
N-body simulations of CDM. By contrast, the TIS model is 
based on a set of physical assumptions concerning the ori- 
gin, evolution, and equilibrium structure of halos in CDM 
universes. Observations of gravitational lenses have the po- 
tential to distinguish between the TIS profile and singular 
ones like the NFW profile, as several observable properties 
of gravitational lenses will be strongly affected by the pres- 
ence, or absence of a central cusp in the density profile. One 
example of an important observable that can distinguish be- 
tween various density profiles is the parity of the number of 
images. Lenses with nonsingular density profiles, such as the 
TIS, obey the odd number theorem. The number of images 
of a given source is always odd, unless the source is extended 
and saddles a caustic (see ISchneider. Ehlers. fc Falcolll992l 
hereafter SEF, p. 172). Lenses with singular profiles, like the 
singular isothermal sphere, the NFW profile, or the Moore 
profile, need not obey this theorem, even for point sources. 
Most observed multiple-image gravitational lenses have ei- 
ther 2 or 4 im ages, and this may argue against profiles with 
a central core jRusin fc Mal200 ll) . There are, however, other 
possible explanations for the absence of a third or fifth im- 
age. That image tends to be very close to the optical axis, 
and might be hidden behind the lens itself. Also, it is usually 
highly demagnified, and might be too faint to be seen. 

We can use the TIS solution to model observed gravita- 
tional lenses individually. Alternatively, we can use the ob- 
servations collectively to constrain the distribution of halo 



properties as characterized by the TIS solution. These prop- 
erties, core radius, velocity dispersion, central density, and 
so on, depend upon the mass of the lensing halos and the 
redshift at which they form. Observational constraints on 
the statistical distribution of these properties will, in turn, 
impose constrains on the cosmological parameters and the 
primordial power spectrum of density fluctuations. 

In this paper, we derive all the lensing properties of 
the TIS. We also compare the TIS with three other density 
profiles: The Navarro-Frenk- White (NFW) density profile, 
the Singular Isothermal Sphere (SIS), and the Schwarzschild 
Lens 1 . To compare the lensing properties of these various 
lens models, we focus on one particular cosmological model, 
the currently favored COBE-normalized ACDM model with 
Q = 0.3, A = 0.7, and H = 70kms _1 Mpc" 1 (this model 
is also cluster-normalized). 

The remainder of this paper is organized as follows. In 
§2, we describe the TIS density profile and the various com- 
parison profiles. In §3, we derive the lens equation. In §4, 
we compute the critical curves and caustics. In §5, we study 
the properties of multiple images: separation, magnification, 
brightness ratios, and time delay. In §6, we study the proper- 
ties of weak lensing, focusing on the average shear. In §7, we 
discuss the likelihood of actually observing cases of strong 
lensing produced by TIS halos. Summary and conclusion are 
presented in §8. 



2 THE DENSITY PROFILES 

In order to compare the predictions for halo lensing for dif- 
ferent halo density profiles, we must relate the parameters 
which define one profile to those which define another. For 
this purpose, we shall assume that all halo profiles contain 
the same total mass within a sphere of the same radius. The 
TIS halo is uniquely specified by the central density po and 
core radius ro. The TIS halo has a well-defined outer radius 
r t at which the mass distribution is truncated, enclosing a 
total mass Mt- There is a unique dimensionless density pro- 
file for the TIS if radius and density are expressed in units 
of ro and po, respectively. Since the central density po is pro- 
portional to p c (Zcoii), the critical density at the epoch of its 
formation, z co n, it is also possible to specify the profile by 
the two parameters, total mass and collapse redshift, (Mt, 
-Zcoii), which is equivalent to specifying the pair (ro, po). It is 
customary to define the total mass of other profiles used to 
model CDM halos as M200, the mass inside a sphere of radius 
r200 with a mean density which is 200 times p c (z) at some 
redshift z. For the sake of direct comparison with the TIS 
profile, we will fix M200 and r2oo for halos of different pro- 
files, which amounts to fixing M t and z co \\ for the TIS halo, 
since M t = 1.167Af 2 oo = 772.6p rj], p = 1.8 x 10 4 p c ( 2co ii), 
and r2oo = 24.2ro if z = z co \\. 

For the NFW profile, there is a third parameter, the 
concentration parameter c, which must be specified in addi- 
tion to the parameters M200 and r 2 oo (or, equivalently, M200 
and z). The value of c is not completely independent of the 
other parameters since there is a statistical expectation that 
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c is correlated with M200 and z. However, for any individual 
halo, c is not known a priori. 

In what follows, we will consider two possibilities for 
comparing our TIS lens with other halo lens models. In the 
first case, the assumptions will be made that the lens redshift 
zl — 2coli for the TIS halo, and the concentration param- 
eter c for the NFW halo of the same mass M200 and r2oo 
at that redshift will be the typical value expected from the 
statistical correlation of c with halo mass and the redshift 
of observation of the halo, z b s = zl- In that case, the ha- 
los are fully specified by the values of M200 and zl- In the 
second case, we can relax the assumption that zl = z c 6a- 
This makes the TIS halo a two-parameter model specified by 
(-M200, Zcoii), and it is assumed that the halo which formed 
at some z co ii ^ zl did not evolve between z co \\ and zl- For 
the comparison NFW halo, we take the same M200 and r2oo 
as the TIS halo, but allow the concentration parameter c 
to take any value, not necessarily the typical one for a halo 
of that mass M200 at redshift zl- In what follows, we fo- 
cus primarily on the first of these cases. This will illustrate 
the relative expectation of the TIS and other profile lenses 
in a sense which is conservative in regard to the possible 
strength of the strong lensing effects predicted by the TIS 
halos relative to the others. In our discussion, §7, we will 
then comment further on the effect of relaxing the assump- 
tion that zl — 2 C oii, as discussed above. 



2.1 The Radial Density Profiles 

The density profile of the TIS is obtained numerically by 
solving a differential equation. However, it is well- fitted by 
the following approximation: 



p(r) = Po 



A 



B 



a 2 + r 2 /r§ b 2 + r 2 jr\ 



(2) 



(Paper I, II) where po is a characteristic density, ro is a 
characteristic radius, and A — 21.38, B = 19.81, a — 3.01, 
6 = 3.82. We will compute the lensing properties of halos 
with this density profile, and compare them with the prop- 
erties derived for three comparison profiles. The first one is 
the Navarro, Frenk, and White (NFW) density profile, 



p(r) 



Pnpw 



( r / r NFw)( r / r NFW + !) 2 



(3) 



where p NPW is a characteristic density and r NPW is a charac- 
teristic radius. We will also consider the Singular Isothermal 
Sphere (SIS) density profile, 



p(r) 



0% 
2nGr 2 



(4) 



where ay is the velocity dispersion and G is the gravitational 
constant. This model might not represent actual halos very 
well, but it is a well-studied profile that has important theo- 
retical value. Finally, for completeness, we will also consider 
the Schwarzschild lens, 



p{r) = M Sch 6\r) 



(5) 



where Ms c h is the lens mass and 5 is the three-dimensional 
delta function. 

To compare the lensing properties of these various 
density profiles, we use essentially the same approach as 



IWright fc Brainerdl ll200ff) . The virial radius r%oo of a halo 
located at redshift z is defined, as usual, as being the ra- 
dius inside which the mean density is equal to 200 times the 
critical density p c (z) = 3H 2 (z)/8ivG at that redshift [where 
H(z) is the Hubble parameter]. The mass M200 inside that 
radius is given by 



Ah 



8007rp c (z)ri 00 



(6) 



When comparing the lensing properties of different density 
profiles, we will consider halos that are located at the same 
redshift and have the same value of T2oo- By definition, these 
halos will also have the same value of M200 ■ By stretching the 
terminology, we will refer to M200 as "the mass of the halo." 
This point needs to be discussed. For the Schwarzschild lens, 
M200 is indeed the mass of the halo. The SIS density profile 
drops as r~ 2 at large r, and the mass therefore diverges un- 
less we introduce a cutoff. The halo mass will then be equal 
to M200 only if the cutoff is chosen to be r2oo- The NFW den- 
sity profile drops as r~ 3 , hence the total mass diverges loga- 
rithmically, and this profile also needs a cutoff. The TIS den- 
sity profile drops asymptotically as r -2 , but the TIS model 
includes a cutoff. This cutoff is located a radius rt fa 1.2r2oo, 
and the mass inside the cutoff is Mt = 1. 168 A/200 • In any 
case, a rigorous definition of the halo mass would require an 
unambiguous determination of the boundary between the 
halo and the background matter (such determination exists 
only for the TIS model), as well as dealing with the fact 
that the assumption of spherical symmetry that enters in 
these models most likely breaks down for real halos at large 
enough radii. Treating M200 as being the actual mass of the 
halo is the best compromise. 

For a given halo mass M200, redshift 2, and cosmological 
background model, all density profiles are fully determined. 
The parameters of the various profiles (po, ro, Pnfwi r NFW> 
ay, and Ms c h) are computed as follows. First, we compute 
f"200 using equation © . The characteristic radius of the TIS 
is then given by 



ro = r 2 oo/C200 = r 2 oo/24.2 



(7) 



(Paper I). The value £200 = 24.2 is actually valid only for 
the Einstein-de Sitter model, but the dependence upon the 
cosmological model is rather weak (Paper II), and for sim- 
plicity we shall use the same value for all models. We now 
integrate equation between r = and r — r2oo, and get 



M200 = 



where 



K- 2 



47TPO?'2ooif200 

C 2 

S200 



A — B — arctan ^22. 4- arctan -^j^- 
C200 a C200 



(8) 



= 2.144 (9) 

(see also lChiba fc Ta kahashi 2002, eq. [10]). By combining 
equations (|SJ and JSf, we can eliminate M200 and r2oo, and 
solve for po. We get 

200p c (z)C 2 2 oo 



Po 



1.8 x 10 p c (a) . 



(10) 



3X200 

The characteristic radius of the NFW profile is given by 

r NFW = r 200/c, (11 
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Figure 1. Concentration parameter c versus halo mass, for NFW 
halos located at redshifts = 0.5 and 1, as labeled. 



where c is the concentration parameter. Unlike the case 
for the TIS, where 7-200 Ao is (nearly) constant, the con- 
centration parameter of the NFW profile is strongly depen- 
dent upon the halo mass, redshift, and cosmological model. 
To determine c for a given halo, we use the formalism of 
lEke. Navarro, fc Steinmet3 (l200ll) . According to this for- 
malism, the concentration parameter in the ACDM model 
considered here varies with halo mass for two illustrative 
values of halo lens redshift zl = 0.5 and 1.0 as shown in 
Figure Once the value of c is known, we can compute the 
parameters of the profile. The characteristic density Pnfw 
is given by 



2Q^p c {z) 

PNFW - 3[ln(l + C )-c/(l + c)] ' (12) 

and the characteristic radius is given by equation 1111 above. 
For the SIS, we integrate equation |gj between r = and 
r = r 2 oo, and get 

2uyr 2 oo 



M200 = 



G 



Combining equations @ and H13^ . we get 



a v 



(13) 



(14) 



Finally, A/g c h is given by M2oo- Figure |2]shows a comparison 
of the various profiles, for halos with the same values of r2oo, 
and M 2 oo- 



JO 




Figure 2. Radial density profiles, for 4 different halos with the 
same values of r200 and M2oo- Solid curve: TIS; dotted curves: 
NFW profiles with concentration parameters c = 2 and 10 (as 
labeled); dashed curve: SIS. 



radius r t w 30ro (the actual value is 29.4 for an Einstein- 
de Sitter universe [Paper I], slightly different for an open 
matter-dominated universe, or a flat universe with a cosmo- 
logical constant [Paper II]). One could always set the limits 
in equation 1151 to ±(r? — £ 2 ) 1//2 . However, it turns out that 
the resulting change in E would be small, for all profiles 
considered. To simplify the algebra, we shall assume that 
equation 1151 remains a good approximation out to r — 00. 
We substitute equations in equation 1151 . For the 

TIS, we get 



i(0 = ""Por-o 



B 



(16) 



[This result was also d erived bv lNataraian. fc Lvnden-Belll 
jl997ft . and lllievl <2000D l. For the NFW profile, we get 



Snfw(C) — 2r NPW p NFW 



1 - 



x < 



x 2 - 1 



1/3, 



■. arg tanh 



1 - x 
1 + x 



x < 1; 



i;(i7) 



2.2 The Projected Surface Density 

The projected surface density is given by 

POO 

E(0= / P(r)dz, 



(15) 



where ^ is the projected distance from the center of the halo, 
and 2 = (r 2 — £ 2 ) 1//2 . In principle, this expression can be 
used only if our expressions for p(r) (eqs. H]-|iD) are vau d 
all the way to r = 00. As we pointed out in the previous 
section, most profiles must be truncated at some finite ra- 
dius rt- For instance, the TIS is truncated at a truncation 



x 2 - 1 



1 - 



Vx 2 



: arctan 



x- 1 
x + 1 



x> 1; 



wher e x = £/)" NFW teaitelmanrj||l99d IWright fc Brainerdl 
l2000l) . For the SIS, we get 



2 

2G£ 



Finally, for the Schwarzschild Lens, we get 
S S ch(0 = Msch5 2 (0- 



(18) 



(19) 
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2.3 The Interior Mass Profile 

For spherically symmetric lenses, one important quantity is 
the interior mass A/(£) inside a cylinder of projected radius 
£ centered around the center of the lens. This quantity is 
given by 

Af (0 = 2ir [ E(£')e'< ■ (20) 
Jo 

We have computed this expression for all density profiles 
considered, using equations llbl - 11911 . We also reexpressed 
the results in units of A/200 with the help of the equations 
given in §2.1. For the TIS, we get 

TV M20O 



Af TIS (£) 



2^200-^200 



-B^Jb 2 +e/r 2 - Aa + Bb] 
For the NFW profile, we get 

A/ N Fw(£) = 47rrNFwPNFW3(£/''NFw) , 

where 

2 



(21) 
(22) 



1+ln 



: argtanh 



1 — x , x 
h In - 

1 + x 2 



x < 1; 



1 



(23) 



: arctan 



x —1 , x 

— + ln 2' 



x > 1: 



I Va? - 1 

^Wright fc Brainerdl l200d : IChiba fc Takahashil 120021) . For 

the SIS, we get 

7rAf 2 oof; 



2r 2 oo 



(24) 



Finally, for the Schwarzschild lens, we get 

A/ S ch(£) = A/200 • (25) 

Figure shows a comparison of Af(£) for various halos with 
the same values of 7-200 and A/200. 



3 THE LENS EQUATION 

Figure 0] illustrates the lensing geometry. The quantities 77 
and £ are the position of the source on the source plane and 
the image on the image plane, respectively, a is the deflec- 
tion angle, and Dl, Ds, and Dls are the angular diameter 
distances between observer and lens, observer and source, 
and lens and source, respectively. The lens equation is 



rj 



-£ - D LS a 



(26) 



[SEF, eq. (2.15b)]. Notice that since the lens is axially sym- 
metric, we can write the quantities 77, £, and a as scalars 
instead of 2-component vectors. We introduce a characteris- 
tic length scale £0, and nondimensionalize the positions and 
deflection angle, as follows: 

D LV 



V 



Co' 



(27) 
(28) 




Figure 3. Interior mass profiles, for 5 different halos with the 
same values of r200 and M2oo- Solid curve: TIS; dotted curves: 
NFW profiles with concentration parameters c = 2 and 10 (as la- 
beled); short-dashed curve: SIS; long-dashed curve: Schwarzschild 
lens. 



observer 














< 

< 




>< 




> 
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Figure 4. The lensing geometry: the dots indicate the location of 
the observer, lensing galaxy, and source. £ and r\ are the positions 
of the image and the source, respectively, and a is the deflection 
angle. The angular diameter distances Dl, Dls, and Dg are also 
indicated. 



D L D LS & 
a = — ^-z • 

Dsto 

The lens equation reduces to 
y = x — a(x) . 



(29) 



(30) 



For spherically symmetric lenses, the deflection angle is 
given by 



a(x) — — I — dx — — I x'k(x')cIx' 

X In S cr it X I 



(31) 



[SEF, eq. (8.3)], where re = E/E cr i t is the convergence, and 
E cr it is the critical surface density, given by 



Ecrit — 



c 2 D s 



(32) 



AnGD L D L s 

where c is the speed of light and Dl, Ds, and Dls are 
the angular diameter distances between observer and lens, 
observer and source, and lens and source, respectively. From 
equations 1201 and 1311 1. we get 

M(tox) 



a(x) = 



7r£oEcrita; 



(33) 
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Hence, the deflection angle is directly related to the interior 
mass. For the TIS, we substitute equation (I2H into equa- 
tion 1331 . and set the characteristic scale £o equal to ro- We 
get 

2irporo 



^crit^ 



[Ay/a 2 + ; 



-B\/b 2 +x 2 - Aa + Bb] . (34) 

This result was also obtained bv lChiba fc Takahashil i200ll) . 
We now introduce the dimensionless central surface density, 
or central convergence, k c , defined by 

^^ = »(^-f), (35) 

and use this expression to eliminate E cr i t in equation 1341 . 
It reduces to 

2abn. 



«Tis(a;) = 



(Ab - Ba)x 



[Ay/a?+: 



-By/b 2 +x 2 - Aa + Bb] 



(36) 

For the NFW profile, we substitute equation 122H into equa- 
tion H33L and set £o = r NFW . We get 



a NFW (zJ 
where 



4n s g(x) 



PNFW r NFW 
Ecrit 



(37) 



(38) 



For the SIS profile, we substitute equation 1241 into equa- 
tion d.'-i-Si . and set 



4-ko v DlDls 



c 2 D s 



GScrit 



We get 



«sis(a;) = 



(39) 



(40) 



(SEF, §8.1.4). Finally, for the Schwarzschild lens, we set 
M(£ x) = M Sch and 



We get 

«sch(^) = - 
x 

(SEF, §8.1.2). 



4GM Sch D L D LS 



c 2 D s 



M Sch 

7rE C rit 



(41) 



(42) 



4 CRITICAL CURVES AND CAUSTICS 
4.1 Solutions 

The determination of the critical curves is quite trivial for 
axially symmetric lenses. Equation 133H can be rewritten as 



m(x) = a(x)x 



(43) 



[SEF, eq. (8.3)], where m(x) = Af(£ :r)Mo£ C rit is the 
dimensionless interior mass. Tangential and radial critical 
curves are defined respectively by 



m(x t ) _ a(x t ) 



xt 



(44) 
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Figure 5. Top panel: Radii of the radial critical circle, x r , tan- 
gential critical circle, xt, and radial caustic, y r , versus k c , for the 
TIS. Bottom panel: xt, Xr, and y r versus re s , for the NFW profile. 



d(m/x) 




da' 


dx 




.dx. 



= 1 . 



(45) 



For the Schwarzschild lens and the SIS, the solutions are 
x t — 1, and there are no real solutions for x r . For the NFW 
and TIS models, we have solved equations I44P and 1451 nu- 
merically for x t and x r . The solutions are plotted in Figure|S] 
as functions of k c and k s . Also plotted is the radial caus- 
tic radius y r , obtained by substituting the value of x r into 
equation I3UH . (The value of y r we obtain is actually neg- 
ative, because the source and image are on opposite sides 
of the lens. The actual radius of the caustic circle, then, is 
the absolute value of y r -) Both x± and \y r \ increase rapidly 
with k c or k s , while the value of x r levels off. Notice that 
there is no relation between k c and k s . The former is the 
central convergence of the TIS, while the latter is the con- 
vergence of the NFW profile at some finite radius and its 
value depends upon the concentration parameter c. Figure||j] 
shows the angular radii of the tangential and radial critical 
circles, 6 r — £r/Di and 8 t = £t/DL, in arc seconds as func- 
tions of the mass of the lens, for lenses located at redshifts 
zl = 0.5 and 1. We assume that the source is located at red- 
shift zs = 3. This is not a a very constraining assumption, 
because the lensing properties vary weakly with the source 
redshift for zs > 1. For instance, for lenses located at red- 
shifts zl ^ 1, the critical surface density E cr i t varies by less 
than 20% when the source is moved from zs = 3 to zs — 5. 

For spherically symmetric lenses, multiple images (and 
thus critical circles) are possible only if the central conver- 
gence re(0) exceeds unity (SEF, p. 236, theorem [e]). For the 
Schwarzschild lens, SIS, and NFW profile, the central con- 
vergence diverges, hence these profiles can always produce 
multiple images. But for the TIS, the central convergence 
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Figure 6. Angular radius of the tangential critical curve (top) and radial critical curve (bottom), versus lens mass, for the TIS (solid 
curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). Results are 
shown for lenses located at redshift ze =0.5 (left panels) and 1 (right panels). 



re(0) = re c is finite, and multiple images can be produced 
only if k c > 1. This explains the sharp low-mass cutoff seen 
in Figure [S] for the TIS (solid curves). 

The value of 9t for the Schwarzschild lens is called the 
Einstein radius 6e- It is often used to estimate the charac- 
teristic scale of image features caused by strong lensing (e.g. 
ring radius, radial location of arc, image separations) and 
to estimate the size of the region within which the mass re- 
sponsible for that strong lensing must be concentrated. Since 
lensing halos are not actually point masses, however, the an- 
gular radius # r i ng of the actual Einstein ring which results 
if the source is located along the line of sight through the 
lens center will usually differ from the Einstein radius 9e, 
assuming that the lens mass distribution is actually capable 
of producing a ring. As we see in Figure 6e significantly 
exceeds Qt for all profiles considered (TIS, NFW profile, and 
SIS) for all masses considered. Hence, a mass estimate based 
on assuming that the scale of image features is of order 9e 
will underestimate the actual mass of the lens, unless the 
lens happens to be a Schwarzschild lens. 

A source located behind the lens will produce multiple 
images if y < y r . The angular cross section for multiple 
imaging is therefore 

2 /„. t.\ 2 



\D S J \ D L 



(46) 



fee 



P E 
b 




In Figure |7| we plot the ratio of the cross sections for the 
NFW and TIS profiles, for sources located at zs = 3. At low 



Mass [M s ] 

Figure 7. Ratio of the cross sections for multiple imaging by the 
TIS and the NFW profile, versus lens mass, for lenses located at 
redshift ze = 0.5 and 1. 



masses, M < 5 x 10 15 M© for z L = 0.5, M < 2 X 10 15 M Q 
for zl = 1.0, the ratios are less than unity, indicating that a 
distribution of lenses described by the NFW profile would be 
more likely to produce cases with multiple images than if the 
same distribution is described by the TIS model. This trend 
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is reversed at higher masses, and a TIS of mass M = 10 16 Mo 
at zl ~ 1.0 is 6 times more likely to produce multiple images 
than a NFW profile of the same mass. 



4.2 Illustrative Example 

Using a simple ray-tracing algorithm, we computed the im- 
age(s) of a circular source of diameter Ay = 1, created by 
a TIS with central convergence k c = 4.015. The results are 
shown in Figure |H| for 8 different locations of the source, 
ranging from y = 8.0 to y — 0.0. For each case, the left panel 
shows the source and the caustic circle (y r — 5.640) on the 
source plane, and the right panel shows the images(s), the 
radial critical circle (x r = 3.334), and the tangential critical 
circle (x t = 9.783) on the image plane. For the cases y = 8.0 
and 6.0, only one image appears. At y = 5.4, the source 
overlaps the caustic, and a second, radially-oriented image 
appears on the radial critical circle. At y = 4.8, the source 
is entirely inside the caustic, and the second image splits in 
two images, located on opposite sides of the radial critical 
circle, forming with the original image a system of 3 aligned 
images. As the source moves toward y = 0, the central image 
moves toward x = and becomes significantly fainter, while 
the other images move toward the tangential critical circle 
and become bright, elongated arcs. At y = 0, the two arcs 
have merged to form an Einstein ring located on top of the 
tangential critical circle, while the central image, very faint, 
is still visible in the center. 

In the case of the TIS and the NFW profile, the radial 
caustics separate regions on the source plane according to 
the number of images a source in those regions produces, 
and when a source moves across such caustic, the number 
of images changes. A curve on the source plane that has 
this property is not always a caustic, although in the case 
of the NFW and TIS profiles, it is. A SIS will produce two 
images if y < 1, and one image if y > 1. However, the circle 
defined by y = 1 is not a caustic. When a source moves 
inside that circle, a second image appears at x — 0, but this 
value of x does not satisfy equation 14511 . and in particular 
the magnification does not diverge at x — 0, as in the case 
of a critical curve, but vanishes instead. Indeed, the total 
magnification varies smoothly and monotonically as a source 
moves across the y — 1 circle. For the Schwarzschild lens, 
there is no radial caustic, and this lens always produces two 
images, no matter where the source is. This happens because 
a(x) diverges as x — > (eq. 1421 ). so no matter how large 
the angular separation between source and lens is, some rays 
will always be deflected toward the observer. However, for 
large separations, the second image is demagnified by an 
enormous factor. 



5 PROPERTIES OF MULTIPLE IMAGES 
5.1 Image Separation 

The locations of the images are computed by solving the 
lens equation 1301 . For the Schwarzschild lens, the solutions 
are 




Figure 9. Plot of a(x) (solid curve) and x — y (dotted lines) 
versus x, for a TIS with k c = 5.005 and 4 particular values of y: 
y = (top line), y = 4, y = y r = 7.515, and y = 12 (bottom line). 
Images are located at values of x corresponding to intersections 
between the lines and the curve. Particular solutions, correspond- 
ing to images located on critical curves, are indicated by arrows. 



(SEF, eq. [8.28]). Hence, there are always two images, with 
separation 



(Az)sch = Vr+4. 

For the SIS, the solutions are 
xi,2 =y±l 



(48) 



(49) 



(SEF, eq. [8.34c]). Again, there are always two images, with 
separation 



(Ax) 



(50) 



Xl,2 = 



y ± y/y 2 + 4 



(47) 



It is an interesting property of the SIS that the separation 
is independent of the source location. For the NFW profile 
and the TIS, the lens equation must be solve numerically. 
In Figure [5] we show the multiple image diagram for the 
particular case of a TIS with k c = 5.005. The solid curves 
shows a(x), while the dotted lines show x — y for particular 
values of y. Each intersection between a line and the curve 
corresponds to one image (see equation 1301 ). If y > y r (bot- 
tom line), the source is outside the caustic circle, and only 
one image appears. If y = y r the source is on the caustic 
circle and a second image appears on the radial critical cir- 
cle, at x = —x r . For y < y r , the source is inside the caustic, 
and the second image splits into two images. Finally, for 
y — (top curve), the central image is located at x = 0, 
and the two outer images are located on the tangential cir- 
cle, at x — ±x t . Actually, these two images merge to form 
an Einstein ring. The slope of the curve a(x) versus x is 
equal to K c at x = 0. It is clear from Figure [5] that if we 
lower k c below 1, any x — y versus x line will intersect the 
curve only once; multiple images cannot occur if k c < 1. 
The corresponding diagram for the NFW profile is qualita- 
tively similar. However, in that case the slope of the curve 
a(x) versus x diverges at x = 0, hence there can always be 3 
images, provided that the source is located inside the radial 
caustic 
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Figure 8. Images of a circular source. Each pair of panels shows the source plane in the left panel, with the caustic, and the image 
plane in the right panel, with the radial (inner) and tangential (outer) critical circles. The position y of the source on the source plane 
is indicated. We used k c = 4.014, and a source of diameter Ay = 1. 



In Figure 1101 we plot the separation between the two 
outer images as a function of the source location. The plot 
only extends to y/y r — 1, since larger values of y only pro- 
duce one image. The solid and dotted curves show the sep- 
arations for the TIS and NFW profile, respectively, with 
various values of k c and k s . The separation is fairly insen- 
sitive to the source location, and stays within ~ 15% of the 
Einstein ring diameter Aa; = 2x t for all values of k c and 
k s considered. For the SIS and Schwarzschild lens, y r is un- 
defined. The dashed line in Figure I1UI shows the constant 
separation for the SIS. For the Schwarzschild lens the sepa- 
ration depends upon y, and therefore cannot be plotted in 
Figure [T01 since y r is undefined. However, in the limit of large 
y, it can be shown that the magnification of the faintest im- 
age drops as 1/y 4 (see SEF, eq. [8.29a]). Hence, in practice, 
the second image will be visible only if y is small, in which 



case equation 1481 shows that the separation is close to 2xt, 
since xt = 1 for that lens. 

Therefore, for all profiles considered, the image separa- 
tion is always of order the Einstein ring diameter, indepen- 
dently of the source location. This is particularly convenient 
for theoretical studie s, when the actual source location can 
be ignored (see, e.g.. iMartel. Premadi. fc Matznerll2002l) . 

5.2 Magnification and Brightness Ratios 

The magnification of an image located at position x on the 
lens plane is given by 



">-=)■ 



(51) 
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Figure 11. Total magnification fitot versus source location rj for lenses of masses 10 12 — 10 16 Mq located at zl = 0.5, for the TIS (solid 
curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



(SEF, eq. [8.17]). We computed the magnification of the im- 
ages produced by halos of masses 10 12 , 10 13 , 10 14 , 10 15 , 
3 x 10 15 , and 10 16 M Q located at redshifts z L = 0.5 and 
zl = 1, for sources located at redshift zs — 3. The results are 
plotted in Figures [TTl415l as functions of the source position. 
Figures [TTI and FHH show the total magnification. The dotted 
curves show the results for a NFW profile. As y (or 77) de- 
creases, the magnification slowly increases, until the source 
reaches the radial caustic y — y r . At that moment, a sec- 
ond image, with infinite magnification (1 — d(m/x)/dx — 
in eq. 1511 ) appears on the radial critical curve (for clarity, 
we truncated those infinite "spikes" in Figs. II lj and 1121 . As 
y keeps decreasing, that second image splits into two im- 
ages, and the total magnification becomes finite again, until 
the source reaches y = 0, and an Einstein ring with infinite 
magnification (1 — m/x 2 = in eq. |51|1 appears on the 
tangential critical curve. Of course, these infinite magnifica- 



tions are not physical, since they can only occur for point 
sources. The total magnification is always larger than unity, 
and always larger when 3 images are present. 

The solid curves in Figures [TTI and [T21 show the results 
for the TIS. At low masses, there is no radial caustic (See 
Fig- Eli an d only one image appears. Because of the presence 
of a fiat density core, the magnification is nearly constant 
if the path of the rays goes near the center of the core. For 
instance, for the case M200 = 1O 12 M0, zl = 0.5 (top left 
panel of Fig. lilt , ro = 7.125 kpc, hence, over the range of 
rj being plotted, we are way inside the core. As the mass 
increases, the magnification increases, until a radial caustic 
forms. This happens at M200 = 1.11 x 1O 15 M for zl = 
0.5, and at M 20 o = 4.7 x 10 15 M Q for z L = 1, according to 
Figure^] At this point, the TIS has the ability to form three 
images, and the results become qualitatively similar to the 
ones for the NFW profile. 
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Figure 12. Total magnification fitot versus source location rj for lenses of masses 10 12 — 10 16 Mq located at zl = 1.0, for the TIS (solid 
curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



The short-dashed and long-dashed curves in Figures fTTl 
and 1121 show the results fot the SIS and the Schwarzschild 
lens, respectively. Because of the absence of radial caustic, 
the magnification always varies smoothly with source po- 
sition, and the only divergence occurs at r\ = 0, where an 
Einstein ring forms. 

For a given mass, the Schwarzschild lens always pro- 
duces the strongest magnification, unless the source is very 
close to the radial caustic of a TIS or a NFW profile, where 
a "spike" of infinite magnification forms. At small masses 
(M200 = 10 12 Mq), the SIS produces a stronger magnifica- 
tion than the NFW profile, for all source positions. However, 
at larger masses, we find a different behavior: If the NFW 
profile produces only one image (that is, we are on the right 
hand side of the dotted spike in Fig. 1111 and 1121 . the SIS 
still produces a larger magnification than the NFW profile. 
But if the NFW profile produces 3 images, then the total 



magnification exceeds the one produced by the SIS. As for 
the TIS, at low masses, where only one image forms, the 
magnification is much smaller than for the other profiles. 
But at large masses, where multiple images can form, the 
magnification becomes comparable to the one for the other 
profiles. 

Figures 1131 1141 and 1151 show the brightness ratios 
I/X1//X3I, and \(j,2/fj,3\, where [li is the magnification of im- 
age i, as a function of source position. By convention, image 
1 is the one between the tangential and radial critical curves, 
image 2 is the one inside the radial critical curve, and image 
3 is the one outside the tangential critical curve (see Fig. |H] 
from left to right, the three images are image 1, 2, and 3). 
The ratio |pi //*3 1 , plotted in Figures [TBI and FPTI always goes 
to unity in the limit 77 — > 0, as images 1 and 3 merge to form 
an Einstein ring. For r\ > 0, image 1 is usually fainter than 
image 3 < 1), unless it is located near the radial 
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Figure 13. Brightness ratio / /j.3 1 versus source location r\ for lenses of masses 10 12 — 10 16 Mq located at zl = 0.5, for the TIS (solid 
curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



critical curve, in which case diverges. For the ranges of 
r\ being considered, The ratio is always very close 

to unity for the Schwarzschild lens, and fairly close to unity 
for the SIS. The TIS and NFW profile do not produce an 
image 1 at large 77, where the source is outside of the radial 
caustic. Notice that in these figures no curves are plotted for 
the TIS at small mass, because only one image forms. 



The ratio | / /J.3 1 , plotted in Figure IT51 for a few cases, 
is much less interesting. Central image 2, which can be pro- 
duced only with the TIS or NFW profile, is usually very 
faint unless it is located near the radial critical curve, in 
which case it might be too close to image 1 to be resolved 
individually. 



5.3 Shear 

The total shear 7(2;) of an image located at position x is 
given by 



7 = 



m(x) 



k(x) 



(52) 



(SEF, eq. [8.15]). For the TIS, we substitute equations 1^ 
and 1431 1 into equation 1521 . and get 



7tis(«) 



abn r 



Ab- Ba 



2A 



\A 2 + : 



2B 



yjb 2 +x 2 -b 



A 



+ 



B 



(53) 



Va 2 + x 2 Vb 2 + x 2 

[this r esult was also obtained bv lNataraian. fc Lvnden-Belll 
( 1997)]. For the NFW profile, the shear is given by 
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Figure 14. Brightness ratio / /j.3 1 versus source location r\ for lenses of masses 10 12 — 10 16 Mq located at zl = 1.0, for the TIS (solid 
curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



f g(x), x^l; 

7nfw(» = k s x < 10 (54) 
[y-41n2, x = l 

l) Wright fc BrainerdEoOfT l. For the SIS, the shear is given 
by 

TsisO) = jo ( 55 ) 

(SEF, page 244). For the Schwarzschild lens, the shear is 
undefined at x — because of the singularity in k (eq. |19)'l. 
Otherwise, 

7sch(2;) = s (56) 

Figures 1161 and 1171 show 71 and 73 (the shear of im- 
ages 1 and 3) versus source position r/ for halos located at 
redhsift zl = 0.5 and 1 with various values of M200, and 



sources located at redshift zs = 3 (notice that 72 is much 
less interesting, since image 2 is usually very faint). When a 
pair of curves appear (two curves with the same line type in 
the same panel) , the upper curve corresponds to 71 and the 
lower one to 73, as the shear is always larger for the image 
located closest to the lens. When only one curve appears, it 
corresponds to 73. 

The Schwarzschild lens always has the largest shear, 
followed by the SIS. The shear for the TIS is negligible at low 
masses, but becomes more important as the mass increases, 
and at Af 2 oo = 10 16 M Q it exceeds the shear for the NFW 
profile and is comparable to the shear for the SIS. 

In the limit 77 — > 0, the shear converges to the value 71, 
73 = 1 for the Schwarzschild lens and 0.5 for the SIS. The 
values for the NFW profile and the TIS vary with mass. 
In the case of the TIS, the shear actually drops to for 
low masses, when only one image forms (see case M200 = 
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Figure 15. Brightness ratio |/j.2 / M3 1 versus source location r\ for lenses of masses 10 12 and 1O 16 M located at z L = 0.5 and 1, for the 
TIS (solid curves) and the NFW profile (dotted curves). 




Figure 10. Separation Aa; between the two outer images, in units 
of the tangential critical radius xt, versus source location y in 
units of the caustic radius y r . The solid curves, from top to bot- 
tom, corresponds to TIS with k c = 10, 5, 2.5, and 1.2, respec- 
tively. The dotted curves, from top to bottom, corresponds to 
NFW profiles with k s = 1.0, 1.0, 0.5, and 0.2, respectively. The 
dashed line corresponds to the SIS. Results for the Schwarzschild 
lens are not plotted. 



10 15 M©, zl = 1-0 in Fig.Hnj. In this circular source 

near rj = produces a nearly-circular image near £ = 0. For 
larger masses, the same source produces three images, with 
two of them (images 1 and 3) located near the tangential 
critical circle, away for the center, where the shear is large. 



The Schwarzschild lens, SIS, and NFW profiles always reach 
this limit at small enough rj, for any value of M2oo- 



5.4 Time Delay 

For axially symmetric lenses, the deflection potential ip(x) 
is defined by 



dip 
dx 



(57) 



For a source at location y producing an image at location x, 
the Fermat potential is defined by 



1 2 

ip,y) = -v) - ^( x ) ' 



(58) 



and the time delay Atij between two images located at Xi 
and Xj, of a source located at y, is given by 



cD l Dls 



(1 + z L )[<j>(xi,y) - <j)(xj,y)] 



(59) 



(SEF, eq. [5.44]). Hence, it is trivial to compute the time 
delay once the deflection potential is known. For the TIS, 
we integrate equation I36H . and get 



'/'tis 0*0 = 
-^4oln 



1abn c 
Ab-Ba 

a + \fo? 



\a\Jo? + x 2 - B^/b 2 + x 2 
b+ ^Jb 2 + x 2 



Bbln 



(60) 
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Figure 16. Shear 71 and 73 of images 1 and 3, versus source location 97 for lenses of masses 10 12 - 10 16 M Q located at z L = 0.5, for 
the TIS (solid curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 
When two curves of the same type appear, the top one is 71, the bottom one 73. 



For the NFW profile, we integrate equation J33> and get 



— In 2 — — 2 argtanh 2 \ / , lacl < 1 : 

2 2 & Vl + W 



1 2 M 2 / \x\ — 1 

- In \- 2 arctan \ / - — : 

2 2 V x +1 



(61) 



\x\ > 1: 



jMeneehetti. Bartelmann. fc Moscardinill2003allrl . For the 
SIS, the deflection potential is given by 



and for the Schwarzschild lens, 



(62) 



tpSchix) = ln\x\ 



(63) 



Notice that for general lenses one normally computes the 
deflection potential ip first, and then differentiates it to get 
the deflection angle a, not the other way around. It is the 
simplicity of spherically symmetric lenses that enables us to 
compute a directly using equation 13111 . 

The time delay between images 1 and 3 (the two out- 
ermost images) is plotted in Figures 1181 and 1191 for the 
same range of masses and source location, lens redshift and 
source redshift as in Figure I11H17I Unlike the magnifica- 
tion and shear, the time delay is not dimensionless, and its 
value varies tremendously over the ranges of M200 and 77 
considered, going from hours to millenia. For the SIS, the 
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Figure 17. Shear 71 and 73 of images 1 and 3, versus source location 97 for lenses of masses 10 12 - 10 16 M Q located at z L = 1.0, for 
the TIS (solid curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 
When two curves of the same type appear, the top one is 71, the bottom one 73. 



time delay varies linearly with r\ (SEF, eq. [8.36]). For the 
Schwarzschild lens, the relation is more complicated (SEF, 
eq. [8.30]), but reduces to linear at small 77. 

For the ranges of M200 and r\ considered, the 
Schwarzschild lens always produces the largest time delay, 
while the other three profiles produce similar delays with 
the NFW profile always producing the smallest delay. At 
low masses, the TIS produces one image, and therefore no 
time delay, but large masses it produces a larger time delay 
than either the NFW profile or the SIS. This indicates that 
there is no simple correspondence between the time delay 
and the central slope of the density profile. 



6 WEAK LENSING 

Weak lensing usually refers to the magnification and distor- 
tion of the image of a background source by a foreground 
lens. Unlike strong lensing, weak lensing normally does not 
produce multiple images of single sources. The detection of 
coherent distortion patterns in the sky has been used to con- 
str ain the mass of clusters. The firs t detections were reported 
bvlTyson. Wenk. fc Vaidesl <199fA : iBonnet. Mellier. fc Fort 
jl994) : iDahle. Maddox fc Lilid dl994l) : iFahlman et a" 
l|l994l) . a nd ISm ail et alj (119951), f ollowe d by many oth- 
ers (see iBartelmann fc Schneider! 1200 ll and references 
therein). More recently, the distortion pattern pro- 
duced by individual galaxies has also been detected 
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Figure 18. Time delay |ti3 | between images 1 and 3 versus source location r} for lenses of masses 10 12 - 1O 16 M located at z L = 0.5, for 
the TIS (solid curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



Dell'Antonio & Tvsonl 11996 


[ lEbbelsl ll98St iHudson et alJ 


19981: iNataraian et alJll99ct 


iFischer et alJhood). 



Following the approach of Wrigh t~fe Brainerdl |_2000), 
we use as measure of the distortion produced by a lens the 
average shear 7 inside a distance £ = r2oo from the lens cen- 
ter. In practice, this quantity would be evaluated by averag- 
ing the shear of all images observed inside T2oo , after having 
eliminated foreground sources. We estimate this quantity by 
integrating the shear over the projected area of the cluster. 
The average shear inside radius x is given by 



70*0 = — 
x 



x'~i(x')dx 



m(x') 



k(x') 



dx . 



(64) 



We eliminate m(x') using equation I4M|I . The two terms in 
the integral can easily be computed using equations 1571 and 
1311 1. respectively. Equation 164H reduces to 



7(a) = -j [#*0 - V'CO)] 



a(x) 



(65) 



For the Schwarzschild lens, ^){x) — In |x|, and therefore 
7 diverges because of the term ij)(Q) in equation I65H . For the 
SIS, a(x) — x/\x\, tp(x) = |as|, and equation I65H reduces to 



7sis(a:) = - 
x 



(66) 
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Figure 19. Time delay [ri3| between images 1 and 3 versus source location r\ for lenses of masses 10 12 - 1O 16 M located at z L = 1.0, for 
the TIS (solid curves), the NFW profile (dotted curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed curves). 



We evaluate this expression at r = r2oo, and get 



7srs0"2oo) = 



4007rp e r 2 oo 

3S cr it 



(67) 



where we used equation (12811 to eliminate £200, then equa- 
tion to eliminate £0, and finally equation 11411 to elimi- 
nate a v . For the TIS, we substitute equations 13611 and I60H 
in equation 1651 . and get 



7tis(z) = 



Aa 



2abn c 



(Ab - Ba)x 2 



x 2 

1 + — 

a 2 



2 In 



1 + yj\ +x 2 /a 2 



~Bb 



1 + 



6 2 



2 In 



1 + ^Jl + x 2 /b 2 



(68) 



We evaluate this expression at r = r2oo, or equivalently x — 
= ?72oo = 24.2. It reduces to 



7Tis( r 2oo) = 408.67 



p c r 200 



(69) 



where we used equation HUH to eliminate po- Both 7 SIS and 
1/3 

7 TIS scale like r^oo , or M 2 g . Their ratio is given by 



( Tsis 



(r = r 20 o) = 1.025, 



(70) 



\7tis 

independent of the redshift zl and mass M200 of the lens 
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NFW profile and the TIS profile, versus mass of the halo, for 
halos located at redshifts zl = 0.5 and 1, as indicated. 



Furthermore, this result is also independent of our assump- 
tions that the background model is ACDM, or the source is 
located at redshift zs = 3. These enter only in the calcu- 
lation of E cr it and p c , which cancel out when we take the 
ratio of equations 1671 and 1691 . For the NFW profile, we 
substitute equations 13711 and 16111 in equation I65H . and get 



4k s 

7npw — — — 



In 2 ? 



, x , 2 / 1 — X 

In - - 4 arg tanh 4 / — — 
2 V l+i 



: arg tanh 



1 - x 
1 + x 



(71) 



To evaluate this expression at r2oo, we simply set x = 
r20o/r NFW — 1/c. 

In Figure |2"U1 we plot the ratio 7nfw/7tis versus mass, 
for lenses located at redshifts Zl = 0.5 and 1, and sources 
located at redshift zs = 3. The ratio 7nfw/7sis can be ob- 
tained by shifting the curves down by an amount log 1.025 = 
0.0107. The ratios are larger at lower redshifts zl- They 
also decrease with increasing mass. This reflects the fact 
that as the mass increases, the concentration parameter c 
of the NFW profile decreases, while the ratio ^oo/^o for 
the TIS remains fixed at 24.2. This figure can be quali- 
tatively compared with th e top left panel of Figure 3 in 
IWright fc Brainerdl ll200Cf) . after the appropriate shifting. 
Notice, however, that (1) the ACDM used by these authors 
is not exactly the same as the one we use, and (2) th e se au - 
thors used the approach of lNavarro. Frenk. fc White! (Il997l) 
to compute the concentration parameter, while we use the 
more recent approach of Eke et al. (2001). At large enough 
masses, the average shear produced by the TIS exceeds the 
one produced by the NFW profile. 

We reach e ss ential ly the same conclusion as 
IWright fc Brainerdl ll2000l) . namely that using the av- 
erage shear to estimate the mass of lensing halos can lead to 
considerable errors if the wrong density profile is assumed. 
At small masses, 7nfw/7tis > 1> an< A therefore the mass 
of a TIS would be underestimated if the lens is incorrectly 
assumed to follow a NFW profile. At high mass, the true 
mass of a TIS would be overestimated. 



7 DISCUSSION 

In the previous two sections, we derived the effects of 
strong and weak lensing, respectively. Here, "strong lens- 
ing" refers to case with multiple images, arcs, or rings, 
while "weak lensing" refers to the magnification and shear 
of single images. We can divide the observed cases of strong 
lensing in two groups. The first group contains the "arc 
second" cases: multiple-image systems with image separa- 
tions of order arc seconds, or rings with radii of that or- 
der fsee iKochanek et al.l Il99ci) . 2 In most cases, the lens 
is a single, massive galaxy, with possibly some additional 
contributio n from the environment in wh ich this galaxy 
is located llTur ner, Qs triker. fc Gotd Il984l : see, however, 
IPremadi fc MarteJl2003h . A classic example is Q0957+561, 
the first gravitational lens to be discovered. The second 
group contains the "arc minute" cases, in which the lens 
is an entire cluster of galaxies. These lenses produce mostly 
giant arcs, with radii in the range 15" — 60" (see Table 1 of 
IWilliams. Navarro, fc Bartelmann||l999l) . The most famous 
case is the cluster CL 0024+1654, which produces multiple 
arcs. 

We showed that for all profiles considered, the image 
separation is weakly dependent on the source location (Fig. 
I10II . when multiple images actually form. If we neglect this 
dependence, the image separation AO « 28t can be read off 
the top panels of Figure^] We see immediately that galaxy- 
size objects cannot produce arc-second separations if they 
are described by the TIS or the NFW profile (under the as- 
sumptions described in §2 that z co n = zl for TIS halos while 
c for NFW halos is the typical value for halos at z h s = zl)- 
The TIS does not produce multiple images unless the mass is 
of order of a cluster mass or above. The NFW profile can pro- 
duce multiple images for any mass, but a separation A9 > 1" 
requires a mass of order 10 13 Mq. Even the SIS needs a mass 
in excess of 1O 12 M0 to produce arc-second separations. At 
the cluster scale, all profiles are capable of producing arc- 
minute separations. In this limit, for a given mass, the sep- 
aration is larger for the Schwarzschild lens, smaller for the 
NFW profile (because of the small concentration parame- 
ter at large mass), and comparable for the SIS and TIS. As 
we indicated in the introduction, the TIS and NFW profile 
are applicable to dwarf galaxies and clusters of galaxies, but 
might not be applicable to ordinary galaxy-scale objects be- 
cause baryonic processes are neglected. Hence, the inability 
of the TIS and NFW profiles to produce arc seconds sepa- 
rations with galaxy-size lenses is not a concern. 

There is very important caveat to this conclusion. So 
far, we have assumed that the redshift zl where the lensing 
halo is located is the same as the redshift z co \\ at which this 
halo first collapsed and virialized. As we pointed out in §2, 
this does not have to be the case, however, since in princi- 
ple the halo could have collapsed at any redshift z co n ^ zl- 
Consider for instance the TIS. The central density po is de- 
termined by equation (111)1 , and so far we have assumed that 
p c (z) was evaluated at zl- But if z co \\ > zl, then p c (z) must 
be evaluated at z C oii, instead. By combining equations ©, 
0, dXHI , and lO, and get 



See |http: / /cfa- www.harvard.edu/castlesl 
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K r OC 



pgr a 
S cr it(zi 



■ oc 



,,1/3 2/3/ s 
M 200Pc (gcoll) 

Ecrit (zl) 



(72) 



Hence, for a given mass, increasing z co ii will result in an 
increase in k c , and lensing will be stronger. Similarly, a 
particular effect (for instance, production of multiple im- 
ages, which requires k c > 1) can be achieved with a smaller 
mass halo if that halo collapsed earlier. For the NFW pro- 
file, the expressions for the characteristic radius and den- 
sity, equations (111! and 1121 . depend on the concentra- 
tion parameter c. The m odels that describe the depen- 
dence of c on Monn. and z <|Navarro. Frenk. fc Whitdll997t 
lEke. Navarro, fc Steinniet3l20C^l) already take into account 
the difference between z and z co \\. Hence, unlike the case for 
the TIS, we are not free to set the redshift z to any other 
value than zl in equation 112H . However, as we explained 
in §2, the value of c we have used so far, which is a single 
function of M200 and zl, is a statistical average over many 
realizations, and a particular halo might have a different 
concentration parameter. 

Let us now focus on the TIS, and perform a simple cal- 
culation to estimate the probability that a TIS can produce 
multiple images. To compute k c , we substitute equations ©, 
Q, and HUH in equation 1351 1. For the particular cosmologi- 
cal model (ACDM) and redshifts (zs — 3, zl = 0.5, 1.0) we 
have considered in this paper, we get 



2/3 



480.75 



Ecrit (zi/) 



Pc(Zc, 



-1 2/3 



Pc(zl) 



Pc(Zcoll) 



Pc(zl) 



2/3 



Mil' * 



M_ 

0.968 , 
289. 



1/3 



z L = 0.5 ; 
z L = 1.0; 



(73) 



where Mis = M200/IO M Q . Strong lensing requires k c > 1. 
If zl = Zcoii, this condition becomes 



M15 > 



(74) 



1.102, z L = 0.5; 
.0.467, z L = 1.0 

(these values are the location of the cutoffs in Fig.|SJ • For the 
cosmological model we consider, a l-cr density fluctuation 
collapsing at redshift z co \\ = 0.5 (1.0) has a mass of about 
M15 = 2 x 10" 3 (4 x 10~ 4 ). Such "typical" objects will 
not be capable of producing multiple images of a source at 
redshift zs = 3, since the resulting value k c = 0.122 (0.095) 
is smaller than unity. This simply indicates that multiple 
images are not produced by the typical objects that collapse 
at Zooli = Zl, which is certainly consistent with the fact that 
fewer than 100 multiple-image systems have been observed. 

Increasing k c above unity would require an object about 
a thousand times more massive than a typical object at 
the same redshift. Objects of this mass are rare but do ex- 
ist. We can make a simple estimate of how atypical such a 
massive object is. Over most of the mass range of cosmo- 
logical interest (from small galaxies to clusters of galaxies) 
the CDM power spectrum can be roughly approximated by 
a power law P(k) oc k n , where k is the wavenumber and 
n w —2. The rms density fluctuation <5 rms is then given by 
S ms w k 3/2 P 1/2 (k) oc k 1/2 . At a given redshift, different val- 
ues of the wavenumber k correspond to different mass scales 
M according to M oc fc~ 3 . The relation between rms den- 
sity fluctuation and mass scale at fixed epoch is therefore 
approximated by 



Increasing the mass by a factor of 1000 therefore reduces cVms 
by a factor of 1000 1,/6 ~ 3. Because of the reduction in <5 rm s, 
a l-cr fluctuation (8 — S rYns ) at this higher mass will no longer 
collapse by the same redshift (it will collapse later), but a 3- 
a fluctuation (S = 35 rms ) will. Such fluctuations are rare, but 
not vanishingly rare. In Gaussian statistics, the probability 
that a randomly located point in space is inside a 3-cr density 
fluctuation (i.e. S ^ 3<5 rms ) is 1/384. Hence, one of every 384 
halos would be capable of producing multiple images (of 
course, whether any halo actually produces multiple images 
depends on the location of the sources) . This was derived by 
assuming z C oii = zl- As equation (I73P shows, the condition 
k c > 1 can be satisfied with a smaller mass it we assume 
that the halo collapsed at an earlier redshift z co \\ > zl- 

A similar calculation could be carried out for the NFW 
profile. In this case, at fixed mass M200 and redshift zl, 
we are free to choose a concentration parameter c that dif- 
fers from the statistical average corresponding to this mass 
and redshift. If we then combine equations dill . 112II . and 
138t . and ignore the weak dependence of the denominator of 
equation I12H on c at large c, we get 



PNFW r NFW c 2 M 

~E 



:rit (Z£,) 



,1/3 2/3/ \ 
200 P c ( Z L) 



Ecrit (zi,) 



(76) 



oc M 



-1/6 



(75) 



Hence, we can increase the effect of lensing, or reduce the 
minimum mass necessary to produce a particular effect, by 
increasing the concentration parameter. 

The simple calculations presented in this discussion 
were for illustration purpose only. In a future paper, we will 
present a detailed calculation of the expected frequency of 
multiple image systems for comparison with the statistics of 
observed lensing. 



8 SUMMARY AND CONCLUSION 

We have derived the lensing properties of cosmological ha- 
los described by the Truncated Isothermal Sphere model. 
The solutions depend on the background cosmological model 
through the critical surface density E cr it, which is a function 
of the cosmological parameters and the source and lens red- 
shifts, and the TIS parameters po and ro, which are functions 
of the mass and collapse redshift of the halo, and the cos- 
mological parameters. By expressing the surface density of 
the halo in units of E cr it and the distances in units of ro, all 
explicit dependences on the cosmological model disappear, 
and the solutions are entirely expressible in terms of two di- 
mensionless parameters, the central convergence k c and the 
scaled position y of the source. We have computed solutions 
for the critical curves and caustics, the image separations, 
the magnification and brightness ratios, the shear, and the 
time delay. The ability of the TIS to produce strong lens- 
ing (multiple images and rings) depends entirely on « c . If 
k c < 1, only one image can form. If k c > 1, either one or 
three images can form, depending on whether the source is 
located outside or inside the radial caustic. When three im- 
ages are produced, the central one is usually very faint, being 
highly demagnified. Degenerate image configurations occur 
when an extended source overlaps a caustic. Two images 
are produced when the source overlaps the radial caustic, 
while an Einstein ring with a central spot is produced when 
the source overlaps the tangential caustic, which is a single 
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point located at y t = 0. These degenerate cases correspond 
to maxima of the total magnification, which diverges as the 
source size goes to zero. When three images are produced, 
the angular separation between the two outermost images 
depends strongly on k c , but only weakly on the source loca- 
tion. 

For comparison, we derived (or extracted from the lit- 
erature) the lensing properties of three comparison models: 
the NFW profile, the singular isothermal sphere, and the 
Schwarzschild lens. Unlike the TIS, all of these profiles have 
a central singularity, which allows them to produce multiple 
images at any mass, provided that the source is sufficiently 
aligned with the lens. In practice, image separations large 
enough to be resolved can be achieved by galactic-mass ob- 
jects only for the Schwarzschild lens, and by supergalactic- 
mass objects for all profiles. 

This paper focused on the intrinsic properties of in- 
dividual lenses described by the TIS model and the com- 
parison models. It provides all the necessary formulae one 
needs to study gravitational lensing by the TIS, and the 
comparison profiles, in specific cosmological models. We will 
present such studies in forthcoming papers. As an illus- 
tration here, we applied the TIS model to the currently- 
favored ACDM universe, to calculate the central conver- 
gence k c expected for TIS halos of different masses and 
collapse epochs. We found that high-redshift sources (e.g. 
zs ~ 3) will be strongly lensed by TIS halos (i.e. k c > 
1) only for cluster-mass halos, assuming that these halos 
formed at the redshift they are observed. As equation 17311 
shows, the halo mass required for strong lensing can be de- 
creased by increasing the formation redshift of the halo. 
An example of a lens for which z co \\ > zl is discussed by 
IShapiro ^^ltev] ||2000j), who showed that the mass - model de- 
rived bv lT^aon^KoTnanski fe DeirAntonTol il99ST) to explain 
their lensing data for cluster CL 0024+1654 at z = 0.39 is 
very well-fit by a TIS halo with p w O.O64/i 2 M pc" 3 and 
ro ~ 20fo _1 kpc. This central density implies that the halo 
collapse redshift is z co n ~ 2.5 (i.e. z co n 3> z l)- 
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